Negative Josephson coupling in the Kondo strong coupling limit 
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We Show that pair hopping through a Kondo singlet give rise to a negative Josephson coupling. 
Thus, our calculation supports the existance of staggered pair correlations in the strong coupling 
limit of a one dimensional Kondo lattice. 



We use a simple strong coupling limit model, of an 
isotropic single channel Kondo impurity, by which we 
show that BCS singlet-pair tunneling through a Kondo 
singlet state results with a negative Josephson coupling 
between BCS superconducting grains. For a lattice of 
such impurities (see figure-1), our result implies a super- 
conducting state with gap nodes due to the modulation, 
or staggering, of the singlet pair correlation function. To 
avoid possible misinterpretations, we stress that the node 
in the pair correlation function due to negative Joseph- 
son coupling is a node in the pair center- of-mass motion. 
It should not be confused with a node in the relative 
pair state (such as p-wave or d-wave), which is usually 
thought to result from the residual repulsive inter- 
actions between electrons on the screened impurity. We 
remark on the relation of our result to the phase diagram 
of the ID Kondo lattice. 

In recent years, several theoretical approaches to the 
Kondo lattice problem have yielded solutions with the in- 
teresting property of having a modulation (or staggering) 
of the pairing order parameter at the Kondo impurities 
Xj\ That is, the singlet pairing correlation function has 
the form, 



X[x 3 -x r ) = (O^dUx') 



(1) 



where xo is the usual correlation of the corresponding 
operator. Coleman et al. Q have use an unusual mean 
field theory, where a coherent three-body bound state of 
the impurity spins with the conduction electrons leads to 
a state with the superconducting order parameter fully 
staggered, i.e., the pair hopping amplitude from site to 
site is negative. Zachar et al. B have used a Bosoniza- 
tion approach to a one dimensional Kondo lattice. They 
found staggered singlet pairing correlations at a special 
exactly solvable " Toulouse point" , which correspond to 
a large value of the Kondo coupling coefficient. Heid 
et al. H have devised a new Ginzburg-Landau theory 
in attempt to explain the superconducting transitions in 
UPt 3 , UPd 2 Al 3 and UNi 2 Al 3 . The distinctive feature 
of this theory || is a superconducting order parameter- 
which resides away from the center of the Brillouin zone. 
Namely, they have finite center-of-mass momentum or 
staggered pairing. 

In all of the above mentioned theoretical approaches, 
the staggered correlation were discussed in conjunction 



with a certain two-channel Kondo H feature leading to 
odd-w (or odd-in-T) pairing ffl,||. For an odd-w pair 
wave function it may be shownthat, the pair transfer en- 
ergy between two odd-w slabs, i.e., Josephson coupling, 
is negative-alternatively. Yet, the reverse is not neces- 
sarily true; having a staggered singlet pairing does not 
entail odd-u> pairing. 

An undercurrent in all of the above approaches is a 
dominant Kondo singlet resonance, and therefore, should 
intuitively be supported by a simple strong coupling ap- 
proximation. Indeed, Coleman et al. || have focussed on 
strong coupling models in which the bare single channel 
Kondo coupling is turned into an effective multi-channel 
coupling due to strong repulsive interactions within the 
conduction electron gas. 

In this paper, we consider a simple strong coupling 
"quantum chemistry" type of model of a single chan- 
nel Kondo impurity. The ground state on the Kondo 
impurity is modeled as a single particle singlet bound- 
state, of one conduction electron with the local impurity 
spin. This approach can be considered as perturbations 
about the infinite Kondo coupling limit (Jk /t — > oo). We 
investigate the problem of Josephson coupling between 
two BCS superconducting grains through an intermedi- 
ate Kondo impurity site in a Kondo singlet state. 

Our main result is that the Josephson coupling is neg- 
ative across the Kondo impurity, even though there is 
no odd-w pairing. Therefore, we demonstrate that stag- 
gered pairing correlations can be a consequence of the 
Kondo singlet state and not only of odd-w pairing. 
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FIG. 1. 



Consider a relatively dilute version of the one dimen- 
sional Kondo lattice model, where the impurity lattice 
constant is much larger than the electron lattice constant 
(see figure- la). We work in the limit Jk ^ t \—V\, 
where t is the hopping matrix element between nearest 
neighbor electron lattice sites, Jk is the Kondo coupling 
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between one local moment and one electron site (note: if 
the Kondo coupling would be to more than one electron 
site then it will become a multi-channel Kondo model), 
and — V < is an attractive interaction between the elec- 
trons (and is of conventional origin, whose details are not 
important for our purpose here). In this strong coupling 
limit, each Kondo impurity spin bounds a single electron 
in a singlet state. The rest of the electrons in each seg- 
ment between impurities are going into a conventional 
BCS state, and thus can be modeled as a BCS grain (see 
figure- lb). 

The effective Josephson Hamiltonian is 

H'" = -2e 2 £ n, [C- 1 ] . . n ! + ]T //,,., (2) 

ij j 



^2 J ij cos ( # 



A 



where 2erij is the charge on grain j, dj is the capaci- 
tance, Aij is line integral of a vector potential from grain 
i to grain j. [rij, e 10 ] — Sije 19 . 

In a Bose liquid, in the absence of a magnetic field, 
the kinetic energy is unfrustrated in the sense that the 
ground state is real and nodeless. It is therefore neces- 
sary that Jij > 0. Since Jij = s i n (g'!_g ) 1S * ne curren t 
across the junction, Jy is a lattice version of the super- 
fluid density N s ( r), (i.e. the proportionality constant 
between the supercurrent and the velocity) 



Ji 



eN s (r) 



J s (r) 
v(r) 



(3) 



Our aim here is to investigate the possibility of having 
negative Josephson coupling, (or negative superfluid den- 
sity) Jij < 0,due to the presence of Kondo impurities in 
the strong coupling limit. 

We investigate the simplest model, in which there are 
only two superconducting grains. In the usual case, in 
which direct single-electron tunneling between grains is 
responsible for Josephson coupling, J 12 is guaranteed to 
be positive in the absence of spin-orbit coupling. How- 
ever, the situation is much richer if we consider the case 
in which the tunneling is indirect, through a localized 
state H . We consider the effect of a Kondo impurity sit- 
uated between the two superconducting grains. Above 
the Kondo temperature, the effect of spin-flip scatter- 
ing [fiof of the tunneling electrons by the impurity spin, 
(which is equivalent to the effect of spin-orbit coupling), 
can produce J\2 < . Typically though, a superconduct- 
ing transition happens well below the Kondo temperature 
(e.g., in heavy fermion materials). 

To investigate the effect of an intermediate coherent 
Kondo state on the Josephson coupling, we examine a 
strong coupling model. In that case, the ground state 
of the system has a single electron bound in a singlet 
state with the impurity spin. We supplement the impu- 
rity Kondo interaction with a residual repulsive on-site 



interaction, U, in order to investigate the effect of such a 
modification, and also in order to keep track of the dou- 
ble occupied impurity states. For simplicity, we neglect 
direct hopping between the two superconducting grains. 
In the context of one dimensional systems, such as the 
Kondo array (see figure), this neglect is justified since 
the direct Josephson coupling between grains is exponen- 
tially small compared with the coupling to the impurity 
(which is situated in the middle between the left and right 
" grains" ) . 

The impurity spectrum is 



Energy 


eigenstates 


E2 = +jJk 


iim,iw>,iTo> = ^(iti) + ]^T» 


E 1= 




E[=0 + Ui 


\m),\wi) 


Eq = —jJk 


|5i> = ^([^>-l^t» 



where ^ is the impurity spin and f is a conduction elec- 
tron spin. The Hamiltonian consists of the two super- 
conducting grains, Hi [a] and H2 [b], the Kondo impurity 
site interactions (Jk and U{), and hopping interaction, 
T, between the impurity and each grain. 

H = H +T (4) 
H a = Hi (a) + H 2 (b) + [j K r • S + Ui (n f] (5) 

T = E [ Tfe ( b l c °* + + h-c.] (6) 



The lowest order contribution to the Josephson cou- 
pling, through intermediate excitations of the impurity 
state, is forth order in the hopping interaction, 



J(4) = j,_ 
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v — Ho + is w — Ho + is w — Ho + is 

JC0S( 0l-0 2 ) 

= - (V>i (a) ^2 J (4) |Vi («) ^2 (&)) \Si) 

= + (Si\ {M2\T±-T±-T±-T \i>ii, 2 ) \Si) . (7) 
-"0 no no 

Where \Si) = ^= (|fr|) - |J||)) is the impurity singlet 
ground state, and the tpj are BCS ground state wave 
functions on grains 1 and 2 respectively, 



1^} = n 



' J i c jq1 c l-ql 



|0> 



(8) 



It is necessary to define a case convention (i.e., the 
order of spin up/down creation operators). We chose to 
put up-spin to the left of down-spin operators in the BCS 
wave- function (||). Thus the proper ordering of the oper- 
ators that transfer a pair from the left grain to the right 

grain is {ijjiip 2 \ (b^b^a-kiak^ \ipiip2>- Direct hoping 
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between grains preserves this order, and results with pos- 
itive Josephson coupling, J > 0. As we shall see, in tun- 
neling processes through the Kondo impurity site, the 
ordering can get permuted, leading to a (— ) sign. This 
effect can be traced to the single electron occupancy of 
the impurity ground state || . 

To fourth order in the hopping interaction, there are 12 
different processes that transfer a pair from the left grain 
to the right grain (and similarly 12 Hermitian conjugate 
processes that transfer from the right to the left). The 
12 processes can be grouped into six pairs of processes 
with identical amplitudes (one of which has the interme- 
diate impurity spin up, fr, and the other process has the 
intermediate impurity spin down, JJ. ). The contribution 
of each and every one of the processes is negative , lead- 
ing to an effective negative Josephson coupling across the 
Kondo impurity. This is our main result. 

rff = (j-l + j 2 + j 3 

= -MY J \T k T q 

kq 



r 4 + J 5 + J 6 ) 
v*u* k v k u q ) 



(9) 
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£k+£q ) \Jit+£ q 



(e k + ^J K ){JK+e k +e q ){lJ K + Ui+e k ) 
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s k +s q 



J K + Ui+e k 



(11) 

(Note that, by definition, all the denominators must be 
positive since we are considering virtual excitations rel- 
ative to the unperturbed ground-state). The first pro- 
cess (see table (1A) in the appendix), is identical to the 
negative Josephson coupling derivation in reference ||, 
which inspired our calculation here. That process would 
dominate in the limit of large repulsion, Ui — > oo, on 
the impurity site (as indeed was the case considered in 
reference In contrast, we consider the limit where 
Ui <C Jk can take any negligible small value. In that 
case, all the 12 processes have contribution of the same 
order of magnitude. In the appendix we list the con- 
tributing processes, and show in two explicit examples 
how the negative sign comes about. 

In conclusion, for antiferromagnetic Kondo interac- 
tions, Jx/t ^ 1) all the contributions to the Josephson 
coupling, from intermediate hopping through the impu- 
rity states, J e M = J1 + J2 + J3 + J4 + J5 + J6, are negative. 
The residual repulsive interactions on the impurity site, 
Ui, have only a quantitative effect. 

In a continuum limit, our result can be associated with 
a recent one dimensional field theoretic derivation by 
Salkola et al. |l2j ; that a local electronic state of charge- 
le spin-0 leads to a ir— kink in the pairing order param- 
eter. 



For a periodic one dimensional Kondo lattice of such 
impurities (figure-1), in the strong coupling limit, the 
pair correlation function is staggered, due to negative 

Josephson coupling across impurity sites, \ (xj — ay ^ = 

(—1)^ J ) xo (xj — ay ^ . Where xo has power law decay 

(as in a ID negative U Hubbard model). 

The form of the staggered BCS pairing, written above, 
is identical to the one which we found previously in a 
Bosonization approach to the ID Kondo lattice P,pT|]. 
That Bosonization calculation was valid in a limited re- 
gion of parameter space. The fact that we find the same 
staggered pairing state in a strong coupling limit (in this 
paper) leads us to propose that a phase of staggered BCS 
pairing correlations extends over a wide region of param- 
eter space in the phase diagram of the ID Kondo lattice 
model. Note that the staggered pairing order, found in 
this paper, is in the strong coupling limit of a single chan- 
nel Kondo lattice, independently of odd — w pairing and 
two-channel physics. 
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A. Appendix: Calculation of Josephson coupling 
terms 

Below, I give a graphical depiction of the 6 processes 
(I A — > 6 A) which transfer a pair from grain — Si to 
grain— S2 , with an intermediate impurity spin— JJ. state. 
(In the initial and final ground states, the impurity site 
is in a singlet bound state; \Si) = -75 (|lU) - l4f)))- Tne 



left column is the intermediate state energy. Since the 
initial and final states are the same for all the processes, 
we indicate them only in the table of the first process: 



(I A) energy 
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(Si) 
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(4A)energy 
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Thus 



' 4 ^ (£k + jJk) (Jk + £k + £ g ) (jJk + U + e k ) 



•J q U k V k Ug) 



I demonstrate the calculation in two cases; In case 
1A, the origin of the (— ) sign is from the operator or- 
dering in the BCS wave function of the grains (i.e., to 

the order of operation (^1^2 1 (^t & l g j. a -fci a fct) l^i^})- 
In case 4A, the (— ) sign results from ordering of the 
electron operators on the impurity site, (1 — tiq s ) riQ a — 

(^ c o&Coa c oa c o<rj ■ ln au of the 12 pair transfer processes, 
one of these fermion operator orderings takes place, lead- 
ing to the negative sign. 
Case-IA: 

Jcos(6»i - 6 2 ) 

= + (s l (ipii>2\ 



E ^ 

kqk q 



-0T a fe'T 
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= +E \ iwi 2 wi UAA^) M) ■ 
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kq 



(Jk + £fe + (| Jk + 



Thus, 



\T k T q \ 2 (v*u* k v k u q ) 



(ek + I Jk) (Jk + £fe + e 9 ) (f Jk 



The factor IfcTq (v*u k VkU q ) is common to all the 
terms. All we need to do is figure out the sign of the per- 
mutation to canonical order of operators and the energy 
denominator of each possible tunneling sequence. 

Case-4A: 



J=+E?l TfcT / (^Tl (0102 
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